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Required Background 



Required Background 

1. Thermodynamics (PHYS253) 

(Covers most of the basic ideas. It would help to look 
through the notes again.) 

2. Quantum and Atomic Physics (PHYS255) 

(You need: Schrodinger's equation, particle in a box, simple 
harmonic oscillator, Zeeman effect, fermions and bosons.) 
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Thermodynamics, briefly. 
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Heat capacities 



This graph shows the heat capacity of copper. At high 
temperatures (T), it tends to a constant. At low T, it falls to 
zero. 




Fig. 10' 2* 4 The molar specific heat of copper: Debye theory with Qd — 309°K 
and experimental points. (After P. Debye, Ann, Physik, vol^ 29, g. 789 (1912).) 

The main theme in Part 1 is to learn how to explain such 
behaviours. 
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Entropy 



The study of heat engines in the 19th century led Rudolf 
Clausius to propose in 1982 the second law of thermodynamics 
in this form: "In a closed cycle, the sum of SQ/T is zero." 

http : //en . wikipedia . org/wiki/Entropy 

SQ is an infinetesimal heat change. Later, SQ/T came to be 
written as SS, or 

dQ = TdS, 

and S became known as the entropy. 

Clausius statement implies that this abstract quantity S is a 
state function. E.g. if you compress and expand a gas and go 
back to the initial p and V at the end, then S also returns to 
the same value. (Note that Q itself is not a state function.) 
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Chemical potential. 



The chemical potential n is equal to the Gibbs free energy at 
equilibrium. The Gibbs free energy is 

G = U + pV — TS. 

One simple use of ji is to find the energy change between two 
phases at constant p and T: 

A/x = AC/ + pAV - TAS. 

At phase transition between liquid and gas, this would be the 
energy difference between 1 mole of liquid changing to gas, and 
1 mole of gas changing to liquid at the same T and p. So 

A/i = 0 

gives the condition for phase equilibrium, because then there is 
no energy advantage in more liquid changing to gas, or the 
other way round. 

http : //en . wikipedia . org/wiki/Gibbs_f ree_energy 
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Quantum mechanics, briefly. 
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Schrodinger's equation. 



This is the ID Schrodinger's equation. 

2m ax 

where m is the mass of a particle, ip is the wavefunction, V{x) 
the potential it experiences, and E its energy. 

In free space, V(x) would be zero. Solving the equation then 
gives 

,0 _ e ikx Qr e ~ikx^ 

where k is related to E given by 

h 2 k 2 

E = . 

2m 
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Momentum. 



The Schrodinger's equation is based on the postulate for the 
momentum p x of a particle: 

-in— = p x . 
ax 

We use this on a wavefunction ip like this: 

—ih— = p x i/j. 
ax 

Applying this to the free space wavefunction e lkx , we find 

p x = hk. 
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Particle in a box 



In a ID box of length L, the potential V(x) could be zero for x 
between 0 and L, and infinity otherwise. 

This means that the wavefunction must be zero outside the 
box, and also zero at the "walls" at x equals 0 and L. 

Inside the box, it is like free space, so ip is e lkx or e ~ lkx . In order 
for ip to be zero at the walls, these can be combined to give 



= sin kx = 



n J/~* rr* n Is* ry* 



2i 

Applying the condition that this is zero at the walls, we find 

k = 



L ' 

where n is an integer. This shows that the energy is quantised 

„ H 2 k 2 n 2 7T 2 H 2 
E = 



2m 2mL 2 
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Statistical mechanics 
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Microstates 



In quantum mechanics, energies are often quantised. We have 
seen that this happens for a particle in a box. 

Consider a system with TV particles. This could be atoms in a 
solid. For now, think of each one like a particle in a box. 

Suppose every particle has energy levels e\, S2, ■■■ 

A microstate refers to a specific arrangement of particles 
among the energy levels, e.g. particle 1 in S3, particle 2 in e\, ... 

A macrostate refers to a distribution of particles among the 
levels, e.g. 6 particles are in e lf 4 particles are in £ 2 > ■■■ 
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Macrostate. 



Suppose there are ni particles in e lf n 2 particles in e 2 » ■■■ 

The macrostate can be written as (ni,n2,...), or simply 
where i = 1, 2, ... Then 

N = ni + 712 H~ ••• 

Suppose the system is thermally isolated from the surrounding, 
so that the total energy E is constant. E is given by 

E = n\E\ + 712^2 "I" ••• 

These two equations impose constraints on the number of 
particles possible on each energy level. 
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Postulate of equal a priori probabilities 



More than one microstate can have the same macrostate. E.g. 
if the macrostate has 4 particles are in e 2 » then these 4 
particles can be any of the TV particles. 

Basic postulate of statistical mechanics: 

"A system is equally likely to be in any microstate that satisfies 
the constraints." 

This means that a system is most likely be be in the 
macrostate with the largest number of microstates. We may 
call this the most probably macrostate. 

Using the method of statistics, it can be shown that if TV is 
very large, then the probability of other macrostates is very 
close to zero. (E.g. when you take the average of more 
readings in an experiement, the error gets smaller.) 
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Constraints 



Assuming the constraints 

N = n\ + n2 + ... 

and 

U = rt\E\ + 712^2 H~ ••• 

it is possible to solve for the most probably macrostate using 
the method of Lagrange multipliers. 

Before that, we need to write down the number of possible 
microstates for a the macrostate. This is just the number of 
ways to arrange TV different objects in some boxes, with n\ 
objects in box (level) 1, objects in box 2, and so on. 

Using the method of combinations, the number of 
arrangements is 

TV" 
n\ !n2l... 

The problem is to find the values of that gives the largest Q. 
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Method of Lagrange multipliers. 



The steps are as follows. 

1. Write down the Lagrange function In Q + AiTV + X 2 U. 
Two new variables, Ai and A2, have to be introduced. 

2. Differentiate with respect to each of ni and Ai, A2. 

3. Set all derivatives to zero. Solve the equations for n^. 
The answer is: 

ni = Aexp(A 2 £;) 

where A is exp(Ai). 

There are two unknowns A and A 2 . In principle, we can 
substitute into the above constraints and solve for these. 
Unfortunately, the equations are too complicated. 
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Boltzmann's postulate 



Around 1872, Ludwig Boltzmann postulated the connection 
between entropy S and the number of macrostates Q: 

S = kg In Q. 

This formula allows us to use thermodynamics, e.g. using 

dQ = TdS, 

to find the answer to the most probable macrostate. Assuming 
Boltzmann's postulate, it can be shown that 

so that the distribution becomes 

rii = Aexp(-Si/k B T). 
This is now called the Boltzmann distribution. 
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Paramagnetic salts 
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Paramagnetic salts. 



Heat capacities tend to increase with temperature, as particles 
can reach higher energy levels. 

At very low temperatures, a paramagnetic salt shows a peak in 
its heat capacity. 




Temperature 77 K 

J.W. Stout, W.B. Hadley, J. Chem. Phys. 40, 55 (1964) 

This unusual behaviour is called the Schottky anomaly. 
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Schottky anomaly 



The Schottky anomaly happens for paramagnetic materials in a 
magnetic field. At very low temperatures, there is little 
contribution from electrons and phonons (lattice vibration). 

Consider a material with ions that have spin 1/2. This is called 
a spin 1/2 paragmagnet. When a magnetic field B is applied, 
the energy level of each ion split into two. 



Energy 




s 2 = juB 



e l = -piB 




As = 2juB 



The energy levels are 



e = — /i.B = —fij^B or + hj^B 
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Lets see if we can "deduce" the Schottky anomally. At 0 K, all 
electrons are at the ground state, so the total energy U is 
—Nii B B, where TV is the number of the ions. 

At high temperatures, the thermal energy the electrons is much 
larger than the difference between the two magnetic energy 
levels. Then the electrons are equally likely to be in —^b^ or 
+H B B, so U is zero 




A sketch of U against T would look like this. 
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Heat capacity 



In principle, we can now find the heat capacity C using 

c = ^. 

dT 

However, the gradient of U neat 0 K is not known. We can find 
this using the Boltzmann distribution. 

Near 0 K, nearly all electrons at at the lower state. The very 
small fraction in the higher state is approximately given by the 
Boltzmann factor exp(—s/k B T) (why?), where 

e = 2/i B B. 

The total energy is then 

U = -Nn B + Nn B exv(-s/k B T). 

Differentiating this shows that the gradient goes to zero as T 
goes to 0 K. 
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Now that we know the gradient is zero at both 0 K and high T 

0 T 




T 

We have just deduced the Schottky anomaly. 
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We shall now use the Boltzmann distribution to derive a 
formula for U and C. There are two levels: 

The total energy is 

The number of electrons in level i is given by 

ui = AexQ(-Si/k B T). 
We must solve for A first. The total number is 

N = A^exp(- e< /fe B T). 
The sum is called the single particle partition function: 

Solving for A gives 

Z SP 
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So 

rii = - — exQ(-ei/k B T). 



The total energy is 



U = J2 n i £ i 



Substituting n % 



U = -z^22 e i exQ(-ei/k B T). 

Z SP 



Substituting for sf 

N 

U = —[-fi B Bexp^ B B/k B T)+LL B Bexp<i-LL B B/k B T)] 

Z SP 

and in Z SP as well: 

_ N -n B Bexp(n B B/k B T) + n B Bexy(-ii B B/k B T) 
exv{ti B B/k B T) + exp(-ti B B/k B T) 
Simplifying, we find U: 

U = -N n B Btanh{n B B /k B T). 
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The formula 

U = -N n B Btanh(iJL B B /k B T). 
does indeed give the graph that we deduced earlier on: 




We note here that a useful formula from the above formula for 
rii, is the probability that an electron is at level sf 

rii exQ(-Si/k B T) 

Pi — — — • 

N Z 
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Differentiating U gives the heat capacity: 

C = Nk B f^'sech 2 (^). 

\k B TJ \k B Tj 

which also gives the Schottky anomaly that we deduced: 
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The entropy. 



We can find the entropy using 



dS = dQ = dU 

T T 

In our case, Q is just U since no work is done. So 

, dU f CdT 

o — 



T J T 

This can be done analytically using the earlier formula for C. 
The result is 



S = Nk B \n 



2 cosh 



Nll b B . ( ll b B 
' -tanh 1 ' B 



T 



\k B T 
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The entropy. 



The graph for the entropy looks like this: 



Nk B ln 2 




Again, we would be interested in the high and low temperature 
behaviours. 



Recall the formula for entropy: S = k B In Q. 
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Behaviour of the entropy 

At low temperature, most atoms would be in the ground state. 

+HbB 

-m,B —tf— 

There is only 1 way to arrange this. So the entropy would tend 
to S = k B \n 1 = 0. 



At high temperature, the difference between the energy levels 
would be small compared with the k B T in exp(—e/k B T). 

+n B B — f — 

-H B B — I 

The atom is equally likely to be in either level. There are 2 
possible arrangements for each of the TV atoms - 2^ in total. 



So the entropy would tend to S = k B In 2 iV = Nk B In 2. 
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Ideal Gas 
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A particle in a 3-D box 



In an ideal gas, we assume that there is little interaction 
between atoms. 

From the kinetic theory of gases, we know that the average 
energy of each atom in an ideal gas is 



1 o 3 

'b 



mv = —kuT. 



2 2 
This means that the total energy is 



U = ^Nk B T 



and the heat capacity is 



c = «L = 3 - NkB 

dT 2 a 

http : //en . wikipedia . org/wiki/Kinetic_theory 



We shall rederive this here using statistical mechanics. The 
ideas will also be used for electrons and phonons. 
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Consider a particle in a 3-D box. 



y 

a 




Solving the Schrodinger's equation, we get the same 
quantisation condition as in the ID box. For the x direction, 



Likewise for the y and the z. 
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Density of states. 



To find the heat capacity of TV particles, we need the energy 

For the spin 1/2 salt, there are only 2 levels to sum. For the 
particle in the box, there are many levels. 

At room temperature, many of these levels could be occupied. 
This is because for an ideal gas, spacing between energy levels 
is much smaller than the average energy of a particle. (Prove 
it.) 

As a result, a plot of against si may look like this. 
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Density of states. 



This suggests that we may approximate the graph to a curve, 
and the sum to an integral. 




Within a small interval de, ni is nearly constant. Then the sum 
of the energies of particles is approximately times the 
number of states in this interval. 

If we divide this number of states by de, we can think of the 
the answer as a kind of density. This density of states (DOS) is 
often denoted by the function g(si). 

The number states in de is then g(ei)de, and the energy in the 
interval given by n^s^(sjds. 
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Density of states. 



The total energy 
can be written as 

Approximating this as an integral gives 

U = J n(s)sg(s)ds. 
where we have replaced the subscript with the argument s. 

This creates a new problem. We must now find the DOS. 
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Counting states 



The wavevectors 



a 

are quantised at uniform intervals in all three directions. 

we can imagine a k space in which the the x coordinate is k x , 
and so on. If we use a point to represent each state, we would 
get a lattice like this. 
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Counting states 




First, we find the number of all states below energy e. This 
means the number of points within 1/8 (why?) of the sphere 
of radius k, where k is the wavevector for e. 

Since the spacing between points is n/a, the volume associated 
with each point is (7r/a) 3 . Therefore, the total number of 
states is: 
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DOS 



We find that the total number of states below k is 

where V is the volume a 3 of the box. 

Suppose k is increased by dk. Then the number of states 
increases by dG. So the density of states is 

dG(k) Vk 2 

We have earlier defined the DOS in terms of energy, g(s). 

Important: To convert the variable to energy, it is not a simple 
matter of substituting the formula 

h 2 k 2 

s = . 

2m 
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Changing variables 



The reason is that for the same interval s, the size of the 
interval dk can be different at different energy because the 
relation 

_ h 2 k 2 
2m 

is not linear. (Prove it.) 

The correct way is to use the method for probability density 
function: 

g £ (s)ds = g k (k)dk. 
http : //en . wikipedia . org/wiki/Probability_density_f unction 

The subscripts are added here to emphasise that g(s) and g{k) 
are different functions. 
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Changing variables 



Rearranging, 

dk 

9e(e) = g k (k) 



de 

We can now substitute the relation 

h 2 k 2 



s = 

2m 

to find g £ (e). The answer is 

, . A-mnV . .1 /9 
^) = ^ r (2m£) 1 / 2 

where we have dropped the subscript again, as is the normal 
practice. (Misleading, but you can tell from the argument.) 

For comparison: 



2ir 2 
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Changing variables 



It is also useful to have the formula for total number states 
in terms of energy. 



Since no e or k intervals are involved, we can simply substitute 

h 2 k 2 

s = . 

2m 

The answer is 
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Ideal gas energy states 



In an ideal gas at room temperature, it is very unlikely for two 
or more atoms to occupy the same energy state. So we can 
assume that each state is occupied by only one atom. This 
simplifies the calculation of the macrostate. 

First, we should justify the assumption. 

At a temperature T, we know from simple kinetic theory of 
ideal gas that the average energy of a gas atom is about 
3k B T/2. 

We need to know how many energy states there are below this 
energy at room temperature. 
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Ideal gas energy states 



We can use the formula we have derived: 

G(e) = ^(2 W£ )^ 
by setting e = 3k B T/2. 

Putting in the numbers for one mole of ideal gas at room 
temperature, we find that the number of states below 3k B T/2 
is about 10 30 . 

1 mole of gas contains about 10 24 atoms. 

This means we have about 10 30 10 24 = 10 6 states per atom 
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Ideal gas energy states 



This means we have about 1 million energy states for every 
atom. 



C 
LU 



So it is extremely unlikely that two atoms would ever occupy 
the same energy state. 

We can use this result to derive the energy distribution of these 
atoms, with the help of the Lagrange multipliers. 
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Indistinguishable particles. 



Consider a interval si around the energy e im Suppose there are 
are atoms and gi states in this interval. 

We need to find the number of microstates. We cannot do it 
the same way as in the paramagnetic salt example. This is 
because the particles in the ideal gas are indistinguishable. 

(Although the ions in the salt are identical, they are 
distingushable - because of their fixed positions in the solid.) 

The problem then is the same as that of arranging identical 
balls and gi identical partitions in a row. 
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Macrostate. 



Use the method of combinations, we find 

(ni + gi)\ 



n^.g^l 



possible arrangements. 



To find the total number of arrangements for all intervals, we 
must multiply together the answer for every interval: 



^ = n 



Since the number of energy states gi is much larger than the 
number of atoms we can simplify this to 



Q = TT- 
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Most probable macrostate 



We can now apply the Lagrange multiplier method. 
The answer is 

ni = Agi exp(-Si/k B T), 
which is the Boltzmann distribution again. 

We can find A using 

N = Y,rii. 

To approximate to an integral as before, recall that the number 
of states Qi in the energy interval is given by g{e)d{e). 

Then ni above becomes 

n(e)de = Ag(s) exp(-s/k B T)ds, 

and 

N = aJ g(e) exp(-Si/k B T)ds. 
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Energy distribution. 



The integral 

Zsp = J g(s)exp(-s/k B T)ds. 

is just continuous energy version of the single particle partition 
function that we have used earlier (for the spin 1/2 salts): 

Z SP = ^TexpC-^/Zc^T). 
So we get the same relation 

N = AZ SP , or A = . 

Z SP 

Substituting the formula for g(s), we find 

f 27rmk B T \ 3 / 2 

Z S P = v [-^^) > 



A = 



N( ^ ^ 3 / 2 



V \27rmkpT 
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Maxwell-Boltzmann distribution 



and 



The energy is just the kinetic energy, e = mv 2 /2. We can 
change to variable v using 

n(v)dv = n(s)ds. 
This gives the velocity distribution 

/ 2 \l/2 / m \ 3 / 2 
n{v) = N i^j f — J exp(-mv 2 /2k B T)v 2 

This is the Maxwell-Boltzmann distribution. 
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At the Columbia University in New York, 



Miller and Kusch measured the speeds of the atoms in 
potassium vapour. They agree very well with the formula 




20 
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Fig. 2. Schematic diagram of the apparatus designed 
to measure velocity distributions. 
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Fig. 4. Typical potassium velocity distributions. The vapor 
pressures in the ovens are given in Table IT. Run 31 was made 
with thick oven slits and runs $1 and 60 with thin slits. 



This was published in Physical Review in 1955. 
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When would it break down? 



We have assumed that: "It is extremely unlikely that two 
atoms would ever occupy the same energy state." 



C 

LU 



Recall that we have arrived at this by looking at the number of 
states below the mean energy 3k B T/2, when T is at room 
temperature. 

If T is very small, there may be far fewer states below the mean 
energy. Then it would be very likely for two atoms to occupy 
the same state, and the Maxwell-Boltzmann distribution would 
not be valid. 

Instead, quantum statistics - such as Fermi-Dirac or 
Bose-Einstein - have to be used. We shall learn about this. 
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Partition function. 



That there is a set of formulae that we can use to find the 
energy from the partition function. 

The energy is given by: 



where the many particle partition function Z N for for 
indistinguishable particles is 




(Prove it.) 



Substituting this, we find 



3Nk B T 
U = — and C = 



dU 




2 

as expected from kinetic theory. 



dT 
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Boltzmann's postulate 



Recall that the Boltzmann distribution comes from 
Boltzmann's postulate 

S = kg In Q. 

The ability to derive the correct heat capacities for 
paramagnetic salts, ideal gases, and other systems are strong 
evidence that the postulate is correct. 

Results such as the heat capacity 




can be directly measured on real gases that are nearly ideal, 
such as Helium. 



This provides a means to obtain the Boltzmann constant, 
which turns out to be related to the gas constant: 

h R 

N A 
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Exercises 
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Some exercises 



Exercise 1 

Minimise the function f(x,y) = x 2 + y 2 , subject to the 
constraint x + y — 2 = 0. 

(i) Do it mentally, or by inspection of the graphs. 

(ii) Do it using Lagrange multiplier. 
Do the results agree? 
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Mentally 




1. The contours of f(x, y) are just circles, e.g. x 2 + y 2 = l, 

x 2 + y 2 = 2, ... The circles get larger as the value increases. 

2. The constraint x + y — 2 = 0 is a straight line with 
intercepts x = 2 or y = 2. 
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Answer: Mentally 




1. By symmetry, f{x,y) is highest at the midpoint between the 
intercepts, where x = y = 1. 

2. Therefore, f(x,y) = l 2 + l 2 = 2. 
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Using Lagrange multiplier 



1. Let g{x, y) = x + y — 2. The constraint is then g(x, y) = 0. 

2. The Lagrange function is 

f(x,y) + \g(x,y) = x +y + A(x + y - 2) (1) 

3. The gradient gives: 

1. Differentiate wrt x: 2x + X = 0. 

2. Differentiate wrt y. 2y + X = 0. 

4. Solving with the constraint equation gives A = — 2 and 
x — y — 1, as before. 
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Some exercises 



Exercise 2 

Use your calculator to work out In 10!. Compare your answer 
with the simple version of Stirling's theorem (AMnTV — AT). How 
big must TV be for the simple version of Stirling's theorem to be 
correct to within (ii) 5% (ii) 1% ? 

Answer 

ln(10!) = 15.11 whereas 10ln(10)-10 = 13.03. There is less 
than 5% difference for N = 24 and less than 1% for B = 91. 
Stirling's approximation is quite accurate even for relatively 
small N. 
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Exercise 3 



Exercise 3 

Consider a 10000 distinguishable particles at room temperature, 
298 K. Suppose that each particle has 2 energy levels, 0.01 eV 
and 0.02 eV. Find the number of the particles in each energy 
level. ( Boltzmann's constant is 1.3807 x 10~ 23 J K _1 .) 
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Exercise 3 



Answer 

Since the number is quite large, we assume that the probability 
at each energy level is given by Boltzmann's distribution 

Aexp(-s//c 5 T). 

At 0.01 eV, Aexp(-s/k B T) = 0.6778A 
At 0.21 eV, Aexp(-s/k B T) = 0.4594A 
The total is 0.6778A + 0.4594A = 10000. So A = 8794. 

Therefore the numbers are: 

at O.OleV, 0.6778 x 8794 = 5961; 
at 0.02eV, 0.4594 x 8794 = 4040. 
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Exercise 4 

Calculate the quantity k B T at room temperature, 298 K. Give 
the answer in eV. ( Boltzmann's constant is 1.3807 x 10 -23 J 
K _1 . Electron charge is 1.6 x 10 -19 C.) 

Answer 

In Joules, k B T = 1.3807 x 10~ 23 x 298 = 4.114 x 10~ 21 J. 

In eV, k B T/e = 4.114 x 10~ 21 /1.6 x 10~ 19 = 0.026 eV. 

Note: Since 1/40 = 0.025 is quite close to this answer, k B T at 
room temperature is often quoted as 1/40 eV. 
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Exercise 5 

Consider a large number of distinguishable particles, at 
temperature T.. Each particle has 4 energy levels 0, k B T, 2k B T 
and 3k B T. Calculate the fraction of particles in each energy 
level. 

Answer 

Since the number is quite large, we assume that the probability 
at each energy level is given by Boltzmann's distribution 

Aexv(-s/k B T). 
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Exercise 5 



First, calculate the following: 

at energy 0, exp(-s/k B T) = exp(O) = 1. 
at energy k B T, exp(—s/k B T) = exp(— 1) = 0.3679. 
at energy 2k B T, exp(-s/k B T) = exp(-2) = 0.1353. 
at energy 3k B T, exp(-s/k B T) = exp(-3) = 0.0498. 

The sum is 1.553. The fractions are, therefore: 

at energy 0, 1/1.553 = 0.6439. 
at energy k B T, 0.3679/1.553 = 0.2369. 
at energy 2k B T, 0.1353/1.553 = 0.0871. 
at energy 3k B T, 0.0498/1.553 = 0.0321. 

Notice that fraction drops to nearly 1/10 after just 2k B T. 

In statistical mechanics, we often want to know if an energy 
level is likely to be populated at some temperature. If the level 
is much higher than k B T, then it is unlikely to be occupied. 
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Exercise 6 

Given that the total energy of a system is 

U = AJ2^ exp(-ei/k B T) 

i 

where A is defined by the normalisation from 

N = Aj2exp(-e l /k B T) 

i 

demonstrate that 

B dT 

where 



Z = 5>xp(-s,//c B T) 
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Answer 



First, note that since Z = J]exp(— si/k B T) 



N = Aj2exp(-s l /k B T) = AZ 



Next, 



<9lnZ ldZ U £j , „ ^ 
8T ZdT Z^k B T 2 K 11 B J 

Multiplying the two equations above, 

<9lnZ A _ , 

== ^2 E^ ex P(-^/^T) 
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Rearranging, 



9 <9lnZ 4 _ 
Nk B T 2 ^- = A^eiexpi-ei/kBT) 



which is U. Therefore 



9 d\nZ 

NkuT = U 
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Exercise 7 

The total energy of a paramagnet is 

U = -Nn B Btanh(iJL B B/k B T) 

Show that the heat capacity is given by 

C-Nk B {-) (exp((?/r) + 1)2 
where 0 = 2/^ B B/k B 

Show that the limiting form at high temperature is 
and at low temperature is 

0\ 2 



C = Nk B (^j exp(-0/T) 
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Answer 



The heat capacity is obtained by differentiating 

U = -Nfj, B Btanh(ij, B B/k B T): 

dT ^ B \ k B T 2 J [ \ k B T 

In terms of 0 = 2\i B Bjk B , 



i -2 



C = Nk 



B 



= Nk B 



— ^ cosh (— 

277 [ \T. 

0 \ 2 



= Nk 



B 



(exp(0/2T) + exp(-6>/2T)) 2 
2 1 exp(0/T) 



.77 (exp(0/2T) + exp(-0/2T)) 2 exp(0/T) 



= Nk B 



exp(0/T) 



77 (exp(0/T) + l) 2 
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At high temperature, exp(#/T) —> 1. 

C = NkJ e -\ 2 6XPWT) 



\TJ (exp(0/T) + l) 2 



2 

(l + l) 2 



4 VT, 

At low temperature, exp(#/T) + 1 — ► exp(#/T) 

0\ 2 exp(0/T) 



J7 (exp(0/T) + l) 2 



6>\ 2 exp(0/T) 



T) exp(6>/T) 2 
= Nk B (£) exp(-0/T) 
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Exercise 8 



Show that the partition function of an ideal gas: 

Z= / ^{2meY' 2 ex^-e/k B T)de 

JO h- 5 

is given by: 

f27rmk B T\ 3 / 2 



Z = V 



[You are given that 



n\{2a) Zn a 



—00 



Basic Statistical Mechanics 



75 



First, move the constant factors outside the integral 

AmirV .1/0 r°° 



Z = 



roc 

(2m) 1 / 2 / s 1 / 2 exp(-s/k B T)ds 
J 0 



h 3 

There is no power of 1/2 in the integration formula that is 
given. So we start with the substitution 



£ = X 2 



Then 

4-mnV ,_ .1/0 r 00 



Z = 



Q (2m) 1 / 2 [°° xexp(-x 2 /k B T).2xdx 
h 6 JO 

(2m) 1 / 2 2 / x 2 exp(-x 2 /k B T)dx 
JO 



h 3 

Next, set n to 1 in the given formula: 

'OO 9 _ 2 2 27T 1 / 2 

x z e a x dx = 



—00 (2a)^a 

The lower of limit of this integral is —00. To make this agree 
with the integral for the partition function Z ... 
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We change the lower limit to 0 and halve the right hand side, 
and the formula becomes: 

•OO 2 _ a 2 x 2 7T 1 / 2 



f 0 (2 a) 2 a 

Now comparing with the partition function again 

A-mnV 



Z = 



(2m) 1 / 2 2 J x 2 exp(-x 2 /k B T)dx 



h 3 

the integrals become the same if we set 



2 1 
a = 



k B T 

Substituting this into the formula gives 



Substituting this into partition function intergral, 
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We get 



Simplifying, 



Z = — — (2m) 1 /2 2 _ - (k B T) 3 / 2 



z = v («) 3 ' 2 



we get the formula for the partition function 
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What is the root-mean-square speed of a helium atom at room 
temperature? 

[Given that the relative atomic mass of helium is 4. 
Atomic mass unit u is 1.6605 x 10 -27 kg. 
Boltzmann's constant k B is 1.3807 x 10~ 23 J K _1 . ] 
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First, we find the mean of v z . This is can be obtained by 



N Jo 

We have multiplied v 2 by the number particles at each energy 
interval de, and divided by the total number N. 

The energy of the atom in an ideal gas is just the kinetic 
energy, 



2 

If we rewrite the above integral for v 2 in terms of the energy, 
we get 




1 



2 



mv 




This integral is just the total energy U, i.e. 



2 _ 



2 



V 



U. 



Nm 
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We have previously obtained the energy U from the partition 
function: 

U = ^Nk B T. 

So the mean square speed is 

- 0 2 3 AT7 _ 3k B T 

v 2 = Nk B T = — — 

Nm 2 ' m 

The root-mean-square speed is then 



v 2 = 



'3k B T 



m 
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We can now calculate the root-mean-square speed. The mass 
of the helium atom is 

Au = 4 x 1.6605 x 10~ 27 = 6.642 x 10~ 27 kg. 
Taking the room temperature as 298 K, the answer is 



v~2 = 



3k B T 



m 



3 x 1.3807 x 10~ 23 x 298 
\ 6.642 x 10- 27 

1363 m/s 
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